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Abstract 

It has recently been shown that relaxation of the rotational energy of hot non-equlibrium 
photofragments (i) slows down significantly with the increase of their initial rotational temper- 
ature and (ii) differs dramatically from the relaxation of the equilibrium rotational energy correla- 
tion function, manifesting thereby breakdown of the linear response description [Science 311, 1907 
(2006)]. We demonstrate that this phenomenon may be caused by the angular momentum de- 
pendence of rotational friction. We have developed the generalized Fokker-Planck equation whose 
rotational friction depends upon angular momentum algebraically. The calculated rotational corre- 
lation functions correspond well to their counterparts obtained via molecular dynamics simulations 
in a broad range of initial non-equilibrium conditions. It is suggested that the angular momentum 
dependence of friction should be taken into account while describing rotational relaxation far from 
equilibrium. 



I. INTRODUCTION 



For several decades, computer simulations, theoretical models, and polarization exper- 
iments have been shaping and deepening our understanding of how molecules reorient in 
liquids and solutions. Most of these studies deal with equilibrium molecular ensembles. 
Much less is known about molecular reorientation under non-equilibrium conditions. 

The study of photodissociation in a condensed phase offers such a possibility. Indeed, 
the incipient photofragments are highly non-equilibrium, and their rotational excitation is 
determined by the excited state potential energy surface of the parent molecule. Tradition- 
ally, polarization photodissociation experiments (both in the steady state and in the time 
domain) have been performed in the gas phase under collision-free conditions.-*^ Recently, 
a number of "real time" measurementa^'^'^'^'^'^'^i'^'^'^'^'^'^'^'^'^'^'^'^'^ and computer 
simulation82^i2ii2^i2^i2i has been reported, where the anisotropy decay of diatomic photo- 
products was studied in the condensed phase. Both classicaF^'^^ and quantun>2S models of 
the anisotropy decay in the dissipative ensemble of photofragments have also been developed. 

So far, theoretical efforts have primarily been focused on studying experimental observ- 
ables, i.e. polarization time-dependent transients (anisotropies). Very recently, an account 
of combined theoretical and experimental study of rotational and orientational relaxation of 
linear CN photofragments under highly non-equilibrium conditions has been published.— 
The authors report a number of extremely interesting and sometimes unexpected phenom- 
ena, which are quite difficult to comprehend within the existing theoretical paradigms. The 
authors demonstrate, in particularly, that relaxation of the rotational energy of hot non- 
equlibrium photofragments (i) slows down significantly with the increase of their initial 
rotational temperature and (ii) differs dramatically from the relaxation of the equilibrium 
rotational energy correlation function, manifesting thereby breakdown of the linear response 
description. What are the physical origins of the long rotational relaxation times and the 
linear response failure? These are the questions which interest us here. We suggest that 
the angular momentum dependence of the rotational friction is responsible for the observed 
behaviour. We offer an explanation why this effect does not manifest itself under equilibrium 
conditions, while it becomes pivotal far from equilibrium. 

The structure of our paper is the following. The analysis of the problem of rotational 
relaxation within the exact generalized master equations, as well as the formal solution of 
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these equations in terms of eigenfunctions is presented in Sec. 2. The generahzed Fokker- 
Planck equation (FPE) with the angular momentum dependent friction is developed in Sec. 
3. In this section, the FPE with algebraic friction is solved for various rotational correlation 
functions (CFs), which are compared with the CFs simulated in Ref.^^ Our main findings 
are summarized in the Conclusion. 

Note that the reduced variables are used throughout the article: time, angular momentum 
and energy are measured in units of a///(A;bT), y/IksT and fc^T, respectively. Here ks is 
the Boltzmann constant, T is the temperature of the equilibrium bath molecules, and I is 
the moment of inertia of the photofragment, so that = ^^/l /{ksT) is the averaged period 
of its free rotation. For CN at 120K, =0.3ps. 

II. GENERAL EQUATIONS 

Let us consider an ensemble of photofragments coupled to a heat bath. We assume that 
the photofragments interact with the bath molecules, but do not interact with each other. 
We suppose that at the time moment t = the joint probability density in the photofrag- 
ment + bath phase space can be written as a product of the equilibrium distribution for the 
bath molecules and a certain non-equilibrium distribution p„e(J) for the photofragments.— 
Restricting our consideration to linear/spherical photofragments, we can apply the projec- 
tion operator technique to the many body Liouville equation and derive the exact rotational 
master equation for photofragmenta^'^'^ 

dtp{3,n,t) = -iAp{3,n,t) - [ dt'C{t-t')p{3,n,t'). (1) 

Jo 

Here p(J,f2,t) is the probability density function, J is the angular momentum of the 
photofragment in its molecular frame, O are the Euler angles which specify orientation 
of the molecular frame with respect to the laboratory one. The free-rotor Liouville operator 
describes the angular momentum driven reorientation, 

A = JL, (2) 

L being the angular momentum operator in the molecular frame. The relaxation operator 
obeys normalization 
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j djdnc{t) = 0, (3) 

which insures the conservation of probabihty (/ dJdQ,p{J, fl, t) = 1) and the detailed balance 

Cit)pBi3) = 0, (4) 

which is responsible for bringing the system under study to the equilibrium rotational Boltz- 
mann distribution 

PB(J) = (27r)-^/2exp{-jV2}. (5) 

Here d = 2 for linear rotors and = 3 for spherical tops. Since molecules are massive inertial 
particles, the relaxation operator C(t) can depend upon L, J, and dj, but (in the absence of 
external fields) is il-independent. This later requirement is a direct consequence of isotropy 
of the rotational phase space. It is important for the further consideration that C{t), by 
its construction, is independent of the initial condition p(J, Q,t = 0). To put it differently: 
once C (t) is chosen, it should describe the time evolution of the probability density function 
p(J, fl,t) for any initial condition. 

We shall further limit ourselves to the long-time (Markovian) evolution of the probability 
density, i.e. assume that the timescale of interest is much longer than the characteristic time 
of bath-induced fluctuations. Thus, the master equation ([1]) transforms into 

dtp{3, n, t) = {-lA - c}p{3, n, t). (6) 

Here 

POO 

C= dtC{t). (7) 
Jo 

In addition, we concentrate on the evolution of the quantities which depend on the angular 
momentum J but are independent of the Euler angles f2. Then, keeping in mind that C 
is independent of il, we can integrate Eq. ([6]) over Q, and arrive at the reduced master 
equation 

dtp{3,t) = -Cp{i,t). (8) 
Here C can depend upon J and 9j. 

Due to the fact that the relaxation operator ([7]) obeys the detailed balance (jlj), 
it is Hermitian. Let us denote its eigenfunctions and eigenvalues by pb(J)\I//j(J) and 
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T]k, respectively {k = 0,1,2,...). The eigenf unctions are assumed to be orthonormal, 
J (iJpB(J)^fc(J)^m(J) = Skm- C possesses a single eigenfunction pB{J)'ifo{J), ^l/o(J) = 1, 
with the eigenvalue rjQ = (this is another way of saying that the Boltzmann distribution 
is the right-hand eigenfunction of C with zero eigenvalue) and all its other eigenfunctions 
PB(J)^fc(J) have positive eigenvalues 1]^. {k = 1,2, ...). If the eigenfunctions are known, we 
can evaluate any CF of interest. Let us assume that, initially, the photofragments had a 
certain non-equilibrium distribution, p„e(J)- Then the angular momentum CF Cj{t), the 
averaged rotational energy Cs(t), and the rotational energy CF CE^t), are determined as 
follows: 

(JJW)ne_g(J^fc(J))ij(J^fc(J)), 

k=l 



CAt) ^ ^ 5^ i™^™^exp{-,,0. (9) 

\ ' ne J,— 1 \ I ne 



\^)ne - \^)b \^)ne " B 

"^'^'^ - {E^)ne'{E)s{E).e " ^ ^^^^^(^M^ ^''^ 

We use the notation (...)^ = / d3pa{3).--, a = B, ne. 

The J-diffusion model^'^'^ and the standard rotational FPEj^i^i^i^ predict the CFs 
Cs{t) and CE{t) to be identical and single-exponential, even in the case of arbitrary initial 
non-equilibrium distribution p„e(J)- The same does the Keilson-Storer model^i^i^, which 
contains the J-diffusion and the FPE models as a special case. On the other hand, the sim- 
ulations carried out in^i^ demonstrate that (i) the averaged rotational energy Cs{t) slows 
down significantly with the increase of the rotational temperature of the photofragments 
and (ii) Cs{t) for hot non-equlibrium photofragments differs dramatically from the equi- 
librium rotational energy CF CE{t), manifesting thereby breakdown of the linear response 
description. 

Eqs. fl9l)- f|TTl) offer a clear explanation of the failure of the standard models of rotational 
relaxation to describe the above phenomena. Within the J-diffusion model2ii2^i2£^ the (stan- 
dard) rotational FPEj^'^'^'^ and the Keilson-Storer model^"^"^, the angular momentum 
CF Cj(t) is described by a single eigenfunction \I^i(J); Cs(t) and CEit) are also described 
by a single, but different, eigenfunction, \I'2(J).^'' Thus, irrespective of the initial condition 
p„e(J), the rate of decay of these CFs is the same. 

If the eigenfunctions which contribute significantly into the coefficients (\E'fe(J))^g, 
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(E^fc(J))„g, {^k{J))B, and (E^fc(J))e in Eqs. ©-([II]) differ from each other, then the 
behavior of Cj(t), Cs(t), and CEit) under equilibrium and non-equilibrium conditions can 
be very different. This is an indication that the standard rotational models must be general- 
ized, allowing for the rate of rotational relaxation to be angular momentum dependent. This 
requirement is easy to understand by using both classical and quantum arguments. Indeed, 
let Tcoii be a characteristic collision time. Such a collision can induce transition between the 
rotational quantum states j and j + I with the frequency ujj+ij provided that^"^ 

(^j+lj^coll = {hi/I)Tcoll < 1. (12) 

The Massey parameter f|T2l) tells us that the molecules possessing relatively small angular 
momentum experience rotationally inelastic collisions, while those possessing high enough 
angular momentum experience j-conserving collisions. This means that the collision rates 
must be j-dependent, and this is explicitly assumed in various "fitting laws" , which are avail- 
able in the literature for describing j-resolved cross-sections.— If we are interested in CPs, 
which are averaged over j, then, under equilibrium conditions, the majority of rotational 
states with non-negligible Boltzmann factors are involved in inelastic collisions. Clearly, if 
the rotationally hot photofragments are produced via dissociation, then the contribution due 
to j-conserving collisions will increase and must be properly accounted for. If the classical 
picture is adopted, we can simply state that the higher is the angular momentum, the more 
collisions are necessary to randomize it. For example, the rate of the angular momentum 
change in a binary collision collision is clearly J-dependent.— It is therefore not surprising 
that Gordon in his classical paper— on rotational relaxation suggested modifications of his 
M-diffusion model toward J-dependent collision frequency. This idea has further been elab- 
orated in papers.— i^'^'^i^i^ However, the effects due to the J-dependent collision rates did 
not receive much attention in the theory of rotational and orientational relaxation, which 
deals primarily with ensemble-averaged rotational and orientational CFs. Indeed, if we wish 
to describe these CFs under equilibrium conditions, then (putting aside non-Markovian ef- 
fects) we can always introduce certain effective (averaged over J) collision frequencies or 
relaxation rates. If we describe rotational relaxation in a wide range of non-equilibrium 
initial conditions, these effective collision frequencies and rates are no longer applicable, and 
their explicit dependence on the angular momentum must be taken into consideration. 
On the basis of detailed simulations, the authors of paper— have proposed the following 
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physical interpretation of the slowing down of rotational relaxation of hot photofargments: 
highly rotationally energetic CN molecules push one (or several) argon atoms out of the 
solvation shell and, after that, rotate more or less freely before the restructuring of the shell 
occurs. This microscopic scenario corroborates entirely with our approach. This is just 
another way of saying that highly rotationally excited molecules experience lower friction 
than their less energetic counterparts. 

For the purposes of the description of CFs fl9l)- f|TT|) under non-equilibrium conditions we 
cannot, unfortunately, use the M-diffusion model with J-dependent collision frequency^ 
and related approaches^'^'^'^'^'^, since all these models assume J-conserving (adiabatic) 
relaxation mechanisms. Neither can we straightforwardly generalize the J-diffusion or the 
Keilson-Storer model by introducing the J-dependent collision frequency, since this proce- 
dure would violate normalization (and therefore conservation) of the probability density (Eq. 
([3])). The next Section is aimed at developing the proper description. 



III. GENERALIZED FOKKER-PLANCK EQUATION 

Without any loss of generality, the relaxation operator ([7j) can be cast into the form of 
the generalized FPE operator: 

C = dji{3 + dj). (13) 

Here the generalized friction ^ is, in general, an operator which depends upon J and dj.—^^^^ 
If the friction is constant (^ = ^ = const) then the standard rotational FPE is recovered, 
which yields the single-exponential CFs^i^i^i^i^ 

Cj{t) = exp{-et}, Csit) = CE{t) = exp{-2et}. (14) 

The concept of friction is fundamental for understanding the rotational dynamics in 
solutions.— i^'^'^i^ In the literature, a distinction has normally been made between the "di- 
electric" friction (which accounts for long-ranged interaction of polar solute molecules with 
a polar solvent) and "mechanical" friction (which is responsible for short-range anisotropic 
interactions). Since argon has been used as a solvent in the simulations carried out in^^, we 
further focus on the "mechanical" friction. As has been explained in the previous Section, 
the drawback off all standard models of rotational relaxation is the following: the relaxation 
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rates are angular momentum independent. To circumvent this problem, we employ the FPE 
f|T3|) with the angular momentum dependent friction (see, e.g.,— »^ and references therein) 
and choose 
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(15) 



1 + 6J2^ ■ 



Here the parameters ^, b > are real, and is an integer.— Such a functional form of the 
friction complies with qualitative considerations presented in Sec. II. Furthermore, the same 
expression has been suggested by Gordon^I as an extension of his M-diffusion model and 
similar formulas are used for describing rotational friction of "active" Brownian particles, 
which convert the energy of the environment into the kinetic energy of their motion.— It 
is necessary to emphasize that the parameters ^, b, and A^, by their construction, must be 
independent of the initial conditions. This means that if Eqs. (fT3|) and (fT5l) describe rota- 
tional relaxation correctly, then they must describe CPs in a wide range of initial conditions 
with fixed ^, b, and A^. 

For simplicity, we shall further restrict ourselves to the one-dimensional case. That is, 
we assume that J is one-dimensional vector (hereafter, the vector notation is therefore 
abandoned), and the equilibrium Boltzmann distribution is given by Eq. ([5]) with d = 1. 
Furthermore, we assume that the non-equilibrium distribution can also be written as a 
Boltzmann distribution, but at a different temperature T„e:— 



We have used Maple 9.5 to numerically solve the FPE f|T3l) with friction f|T5l) and calculate 
Cj{t), Cs(t), and CE{t). The results of the calculations are depicted in Fig. 1. As expected, 
the decay times of all the CFs increase when we approach non-equilibrium conditions (hot 
photofragments). All the non-equilibrium CFs Cj(t), Csit), and CE(t) differ considerably 
from each other as well as from their equilibrium counterparts (IT^ . As is clearly seen, 
Cs{t) calculated for hot photofragments deviates significantly from the equilibrium rota- 
tional energy CF CE(t) (compare the upper dashed line with the lower dotted line). This is 
in the agreement with the results of papers,—"^ which demonstrate inadequacy of the linear 
response theory in reproducing Cs{t) far from equilibrium. 

In order to get more insight into the behavior of CFs, it is helpful to consider their integral 
correlation times 



pUJ) = (27r/x)-^/'exp{-xJV2}, X = T/T^,. 



(16) 
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Ta = / dtCa{t), a = J,E, S. (17) 

Now one more argument in favor of choosing algebraic friction ( JTSll is coming: as is shown 
in the Appendix, we can calculate Eqs. (fT7|) analytically. The results read: 



(IS, 



1 ,(2iV+l)!! i/x^+i-n 

2?\ A'+l 3/x-l J ^ ' 

If 6 = 0, then we recover the standard FPE formulas. The presence of 6 7^ causes the 

increase of the relaxation times. 

Eqs. fll8l)- fl20l) help us to reveal an important fact. Let us assume that the J-dependence 

of friction (fT5|) is weak, that is 6 <^ 1. If the initial conditions are close to equilibrium (x ~ 1), 

then the contribution due to the J-dependent friction in Eqs. f|T8l) - fl20l) is proportional to 

b and is therefore small. This observation supports the use of the models with constant 

relaxation rates (the J-diffusion model,— i^*^ the standard rotational FPE,— '^'^'^ and the 

Keilson-Storer model^*^"^) for the description of molecular reorientation under equilibrium 

conditions. Suppose that the initial conditions are highly non-equilibrium, so that the 

characteristic temperature of the photofragments is much higher than the bath temperature, 

X = T/Tne ^ 1- Then the contribution of the second term in Eqs. (fT8|) - (l20|) is determined 

by the factor of h/x^ ■, which is no longer small even if 6 <^ 1.— Thus 

T„e r^T/Vh (21) 

delivers the threshold value of the non-equilibrium temperature (or rotational energy), for 
which the slowing down of rotational relaxation becomes significant. The existence of such a 
threshold temperature is clearly demonstrated in Ref.— The farther we are from equilibrium 
(the higher is x)i the stronger is the effect due to the J-dependence of friction. This finding 
corroborates a general paradigm of statistical physics and non-equilibrium thermodynamics, 
which emphasizes the growing role of fluctuations far from equilibrium. 

The above considerations explain why the linear-response theory, which identifies the 
non-equilibrium averaged rotational energy Cs{t) with equilibrium rotational energy fluc- 
tuations CE{t), breaks down far from equilibrium. If we are close to equilibrium, then the 



contribution of the high-J states into rotational relaxation is minor. Thus the equilibrium 
fluctuations sample adequately the angular momentum exchange between the solute and 
solvent molecules. If we move far from equilibrium, then rotational relaxation spreads over 
the high-J states, and equilibrium fluctuations are unable to adequately sample the angular 
momentum exchange. 

Once Eqs. (fT8|) - (!20|) are available, we can attempt to make more quantitative comparison 
of the present theory with the results of simulations reported in.^^ If our approach grasps 
essential physics of the phenomenon, we shall be able to fit all the simulated Cs{t) by a 
single set of the parameters N, ^, and b. We have thus proceeded as follows. We have 
numerically calculated the integral relaxation times of Cs{t) simulated at different T^e- By 
combining the obtained dimensionless values of ts at equilibrium [T^e = T = 120K) and 
at a certain temperature Tne > T we have used Eq. f|T9|) to calculate the corresponding 
parameters ^ and b for = 1, 2. The results of these calculations are summed up in Table 
1. If our approach is self-consistent, the values of ^ and b must be more or less the same 
for any T„e. As is seen from the Table, this is certainly not the case for = 1, while the 
calculations with N = 2 exhibit a remarkable consistency. Thus we can speculate that the 
FPE f fT3l) with friction f|T5l) at = 2 describes the simulated Cs{t) rather adequately. The 
so calculated Cs(t), along with their simulated counterparts^ are depicted in Fig. 2. The 
overall agreement is satisfactory, but not perfect. There exists a number of reasons for such 
a deviation. First, the simulated CFs clearly exhibit short-time non-Markovian behavior. In 
order to describe this, the present FPE should be generalized to account for memory effects 
(see, e.g.,—). Such a generalization is necessary, for example, for reproducing the short-time 
(~ lOOfs) coincidence of the simulated Cs(t) at different Tne-— Second, in order to calculate 
the integral relaxation times f|T8l) - fl20l) analytically, we have switched to the one-dimensional 
model, while the actual description for linear photofragments must be two-dimensional. 
Third, as is well known, the FPE itself may not be very good for reproducing rotational 
relaxation in liquids^, so that more refine approaches might be necessary. The present paper 
undertakes the first step toward understanding non-equilibrium rotational relaxation. The 
message is as follows: the friction, or the collision frequency, or the relaxation rate must 
be taken J-dependent in order to adequately describe rotational relaxation under highly 
non-equilibrium conditions. 



10 



IV. CONCLUSION 



Recently, the groups of Bradforth and Stratt have pubhshed the results of the combined 
theoretical and experimental study of rotational relaxation of CN-photofragments in a con- 
densed phase.— >^ They show via molecular dynamics simulations that the time evolution of 
the averaged rotational energy, Cs{t), (i) slows down dramatically with the increase of the 
rotational temperature of the incipient photofragments and (ii) deviates significantly from 
the equilibrium rotational energy CF, CE{t), manifesting the linear response breakdown. To 
comprehend and explain this unusual behavior, we have developed a theory of rotational 
relaxation under non-equilibrium conditions, by extending the standard FPE approach to 
account for angular momentum dependent friction. We have calculated the angular mo- 
mentum CF Cj(t), as well as Cs(t) and C_E(t), and compared them with their simulated 
counterparts.™ We have shown that J-dependence of rotational friction is responsible for the 
observed slowing down of rotational relaxation and failure of the linear response description, 
provided the effective non-equilibrium rotational temperature exceeds the threshold value 
fl^T]) . The farther the system is from the equilibrium, the stronger is the effect. 
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APPENDIX A: CALCULATION OF THE INTEGRAL RELAXATION TIMES 

We start from the one-dimensional FPE with the angular momentum dependent friction 

dtp{J,t) = dji{J){J + dj)p{J,t). (Al) 
Evidently, a CF of the functions A{J) and B{J) can be evaluated via Eq. flAip as follows: 

/oo 
dJB{J)p{J,t) (A2) 
-c» 
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(the notation is identical to that used in Eqs. (pi)- (fTT!) ). The initial condition to Eq. (lAip 
reads then: 

p{J,Q) = A{J)pUJ)-{A{J))^^^Ps{J). 

Here the initial non-equilibrium distribution is determined by Eq. fll6p . and the Boltzmann 
equilibrium distribution is defined via Eq. ^ with d = 1. Let 



oo 



TAB = / dt{{A{J)B{jmne- W)))neW))B] 







be the integral relaxation time for the above CP. Then, integrating Eq. (lAip over time and 
introducing the quantity 

poo 

r]{J)= / dtp{J,t), 







we can express tab through 77 (J) as follows: 

/oo 
dJB{J)r]{J). (A3) 
-00 

Here 1]{J) obeys the differential equation 

- {A(J)p„e(J) - {AiJ))^,^PBiJ)} = dji{J){J + dMJ). (A4) 

Eq. ( 1A4[) must be solved with the boundary conditions rji^oo) = 0. It is convenient to 
introduce the quantity rj{J) via the expression 

viJ)^viJ)PBiJ). (A5) 
Upon the insertion of the above equation into Eq. flA4p . we get: 

^i-— f dJ' {A{J')pUJ') - (A(J))„,p^(J')} = d.j?j{J) (A6) 
Pb\J)c,[J) Jo 

(the integration limits have been chosen to comply with the boundary conditions ri{±oo) = 
0). As is seen, Eq. flA6p can be integrated, in principle, for any Pne{J), ^{J), and A{J). To 
simplify the subsequent presentation, we shall limit ourselves to the case when 

B{J) = J"',m = loT2. (A7) 

Inserting this expression into Eq. flA3p . using the definition flA5l) and integrating by parts 
yields then 

POO 

TAB= / dJJ"'-^pBiJ)djviJ), m = lor2. (A8) 



m—l 

UJ J 

00 
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Combining Eqs. (1A6|1 and (lASp . we get: 

f OO -1 /"J 



Tab 



/ dJJ^-'—- dJ'{A{J')pne{J')-{A{J))^^pB{J')}, m = loT2. (A9) 

J-oo SWJ Jo 



Using the explicit form f|T5l) of the angular momentum friction and integrating Eq. flA9l) by 
parts gives, finally: 

POO / jm jm+2N \ 

Tab = t dJ [ Vh 



. ^^^^^j{A{J)p^,{J)-{A{J))^^PB{J)}, m = loT2. (AlO) 
This is the formula which has been used for the calculation of integral relaxation times 

m-m. 
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TABLE I: The dimensionless parameters and b calculated via Eq. p9p from the integral re- 
laxation times Ts (given in units of r^) simulated at different temperatures Tne (in K). Integral 
relaxation time ts = 1.77 at equilibrium temperature 120K. 



Simulation^^ 


= 1 


N = 2 


Tne Ts 


e b 


e b 


1917 5.29 


0.39 0.12 


0.29 0.0015 


2395 7.73 


0.44 0.18 


0.29 0.0017 


2875 11.25 


0.53 0.29 


0.29 0.0018 
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time, ps 

FIG. 1: Angular momentum CF Cj{t) (full lines), averaged energy Cs{t) (dashed lines) and 
rotational energy CF CE{t) (dotted lines) for algebraic friction ()15p with = 1, ^ = 1, and 
b = 0.29. The lower CFs are calculated for equilibrium conditions (x = T/Tne = 1) and the upper 
CFs are computed for hot photofragments (x = 0.1). 
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FIG. 2: Averaged energy CF Cs{t) for Tm = 120K (a), 1917K (b), 2395K (c) and 2875K (d). Full 
lines correspond to the molecular dynamics simulations,— and dashed lines display the results of 
the FPE calculations with algebraic friction (jl5p . N = 2 and the parameters and b are listed in 
Table 1. 
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